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Classical properties of 1/4 BPS dyons were previously well understood both in the field 
theory context and in the string theory context. Its quantum properties, however, have 
been more difficult to probe, although the elementary information of the supermultiplet 
structures is known from a perturbative construction. Recently, a low energy effective 
theory of monopoles was constructed and argued to contain these dyons as quantum 
bound states. In this paper, we find these dyonic bound states explicitly in the N = 
4 supersymmetric low energy effective theory. After identifying the correct angular 
momentum operators, we motivate an anti-self-dual ansatz for all BPS bound states. 
The wavefunctions are found explicitly, whose spin contents and degeneracies match 
exactly the expected results. 
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1 Introduction 



In the N=4 supersymmetric Yang- Mills theories, there can be 1/2 BPS and 1/4 BPS configurations. 
The precise nature of these states depends on the asymptotic values of six Higgs fields in the 
theory. When the Higgs expectation values have only one independent component, only 1/2 BPS 
configurations can appear. Classically, 1/2 BPS configurations are made of monopoles or dyons, 
and its electric field is proportional to its magnetic field. When the Higgs expectation values have 
two or more independent components, 1/4 BPS configurations can also appear, which are all dyons 
whose electric charge is not proportional to its magnetic charge (l|, ||, ||. 

Any 1/2 BPS configuration of a given collection of monopoles is specified by their moduli 
parameters, and the low energy dynamics of these 1 /2 BPS monopoles is determined by the metric 
on the manifold spanned by these moduli parameters On the other hand, as shown in Ref. [^, 
solutions to 1/4 BPS equations can be obtained in two steps: First, one solves for purely magnetic 
soliton, which may be regarded as 1 /2 BPS in a technical sense. Then, one solves for certain linear 
combination of gauge zero modes in this purely magnetic background. Curiously enough, one can 
build the electric part of the 1/4 BPS dyons from such gauge zero mode. Because existence of 
gauge zero modes are guaranteed for any 1/2 BPS monopole [||, any 1/4 BPS configuration is 
again specified by the moduli parameters of the corresponding 1/2 BPS monopole configuration. 
Given a fixed set of moduli parameters, the solution of the second BPS equation determines relative 
part of the electric charge of monopoles uniquely. 

One may regard 1/4 BPS configurations as deformed 1/2 BPS configurations when the ad- 
ditional and independent Higgs expectation is turned on. When this second Higgs expectation is 
quite small compared to the first, the deviation of the 1/4 BPS configurations from the 1/2 BPS 
configuration is small. In such cases, the authors (with C. Lee) have shown in a recent paper j|] 
that one can describe the low energy dynamics of both 1/2 BPS and 1/4 BPS configurations with 
an effective nonrelativistic Lagrangian. 

The kinetic part of the Lagrangian is given by the moduli space metric of the 1/2 BPS 
configurations. The potential is also present, and is given by the square of the norm of a tri- 
holomorphic Killing vector field related to an unbroken U(l) gauge symmetry. The size of this 
attractive potential is proportional to the square of the additional Higgs expectation value. This 
effective Lagrangian can be interpreted as low energy dynamics of 1 /2 BPS monopoles with attrac- 
tive potential, in other words, and the 1/4 BPS configurations should be realized as BPS bound 
states of monopoles with additional electric quantum numbers. 



1 



In Ref. ||, the full N=4 super symmetric low energy effective Lagrangian is written. This is 
a sigma model with potential that has extended complex supersymmetry with a central term, and, 
as usual, the wavefunctions can be interpreted as differential forms on the moduli space. The BPS 
equation was found and translated to the language of differential forms. 

The simplest nontrivial 1/4 BPS configurations appear as composite of two distinct funda- 
mental magnetic monopoles in SU(3) gauge theory. Furthermore, the 8-dimensional moduli space 
of these two monopoles is known exactly. From Ref. ||], several facts are known about this case. 
First of all, classical 1/4 BPS configurations are made of two 1/2 BPS dyons at rest, whose mutual 
distance is determined by their relative electric charge. Also, the supermultiplet structures of all 
such dyons have been found by the perturbative method around the 1/2 BPS state of the zero 
relative electric charge, where the nonzero relative charge states are constructed by exciting certain 
massive excitations on 1/2 BPS configurations. 

In this paper, we reconstruct these 1/4 BPS dyons as quantum bound state of two distinct 
SU (3) monopoles in the low energy dynamics described above. We construct all such SU (3) dyons. 
We also recover the phenomenon of instability found in Ref. []l], [|]: The bound state wavefunction 
loses its normalizability exactly at the point where the instability should set in. Furthermore, we 
explicitly show that each (stable) dyon comes in the same supermultiplet as found in Ref. Q 

The plan of the paper is as follows. In Sec. 2, we briefly discuss the moduli space of a pair of 
distinct monopoles in SU(3) theory. In Sec. 3, we review briefly the supersymmetric Hamiltonian 
and BPS conditions on wavefunctions as shown in Ref. ||. In Sec. 4, we discuss the angular 
momentum and an ansatz for the BPS wave functions. In Sec. 5, we solve the BPS equations. In 
Sec. 6, we conclude with some comments. 

2 A Pair of Distinct Monopoles in the SU(3) Gauge Theory 

Consider N = 4 SU(3) gauge theory spontaneously broken to U(l) 2 . When the six Higgs expecta- 
tions are all collinear, the theory contains two distinct types of fundamental monopoles, which we 
will label by a. and j3. The low energy interaction between a. and (3 monopoles can be described 
by the moduli space dynamics. There are four collective coordinates for each monopole, three for 
its position and one for the U(l) phase. We call their positions and phases to be Xi,Xi,i = 1,2, 
for ot and (3 monopoles, respectively. Let us parameterize the masses of these monopoles as n\ and 
H2- We are suppressing the gauge coupling constant in all subsequent formulae. 
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The exact nonrelativistic effective Lagrangian has been found to be a sum of the Lagrangians 
for the center of mass and the relative motion Q. As there is no external force, the center of mass 
Lagrangian is a free one: 

where the center of mass position is X = (/Uixi + ^2X2) /(fii + ^2) and the center of mass phase 
is XT = Xi + X2- The relative motion between them is more complicated and described by the 
Taub-NUT metric || g], and has the Lagrangian, 

where the relative position is r = X2 — xi, the relative phase is \ = 2(^1X2 — / (/^l + ^2), arid 
w(r) is the Dirac potential such that V x w(r) = — r/r 3 . The range of x 1S [0> 47r]. From now on, 
we will suppress the scale \x by setting fj, = 1. The resulting monopole moduli space metric is then, 

9rel =( 1 + l) * 2 + (iTTaO ( ^ + W(r) " * )2 ' (3) 

up to an overall scale. This Taub-NUT space, Aio, has the topology of R A = R + x S 3 . The 
eight-dimensional, total moduli space is then given by 

(4) 

where Z is the identification map 

(XT,X) = (xT + 27r,x + ^^)- (5) 

For later convenience, we will make another choice of coordinates involving Euler angles on 

S 3 , 

ftel = (l + ;) [dr 2 + r 2 a\ + rV 2 2 ] + a 2 , (6) 

where the <7 a 's are 1-form frames on S s , and satisfy the canonical relationship, 

da a = 2 e abc<Tb A a c . (7) 
More explicitly, we may write these 1-forms in terms of SU (2) Euler angles as follows 



o\ = — sin xdO + cos x sin ( 
(J2 = cos xdO + sin x sin I 
^3 = dx + cos (8) 
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The ranges of 9, (j), x are respectively tt, 2ir, Air. Let us define an orthonormal basis uj^ by 



uj° = Jl + l/rdr, 



u 1 = Vr + T(7i, 



uJ 2 = \/ r 2 + r (72, 



- 3 = ^<*- (9) 
Because the Taub-NUT manifold is a hyperKahler 4-manifold, its curvature is anti-self-dual with 
an appropriate choice of orientation. 

When the Higgs vacua is slightly misaligned, the two monopoles are attracted to each other 
10 1 . The effective low energy potential U of this static force has been found in Ref. || for all 



multimonopole configurations in all N = 4 gauge theories. Specializing to the case of a pair of 
distinct monopoles in SU(3), the relative part of this potential is given by a squared norm of the 
Killing vector field d x , up to an overall factor, 

^ rel = Ia 2 /A A\ (10) 



2 \dx dx, 

where a is a measure of Higgs misalignment. The interacting part of the two monopole dynamics 
is dictated by an effective Lagrangian, whose bosonic part is 

Aei = ~ - W reI = \ (l + I) r 2 + \ (j-^) (X + w(r) -rf-\ . (11) 

Note that the potential U ve \ increases from zero at origin to a 2 /2 at infinity. This behavior allows 
new bound states of the dynamics which would not have been possible for a = 0. Among these are 
certain dyonic states that preserve 1/4 of field theory supersymmetries. The purpose of this note 
is to reconstruct these 1/4 BPS dyons as BPS quantum bound states in the low energy dynamics 
of monopoles. 



3 Super symmetry and BPS Bound 

We begin by recapitulating generic properties of the N=4 supersymmetric quantum extension of 
the bosonic effective action || [Uj 12]. Its form is rather similar to the usual supersymmetric sigma 
model action but supplemented by an attractive bosonic potential together with its fermionic 
counter part. These potentials are determined by a single Killing vector field G. The supersym- 
metric Lagrangian written with real fermions is 

-g^ v G^G u - D^^r) , (12) 
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where ip 1 * is a two-component anticommuting Majorana spinor and 7 = 02, 75 = 03, and ifi = ip T ^/°. 
In case of relative dynamics of the two SU(3) monopoles, G is equal to ad x . As required for 
the N=4 supersymmetry, the metric here is hyperKahler, endowed with three complex structures 
iWn v ( a _ 1; 2, 3) that satisfy 

j(a) T (6) = _ S ab + e abc z {c)^ ^ 

D^% = 0. (14) 

For the sake of N=4 supersymmetry, the Killing vector G = a- K should be triholomorphic; namely 
its action preserves the three complex structures via 

C G l {a) = 0, (15) 

where C denotes the Lie derivative. 

Upon quantization, the spinors tfj A = e^ip^ with vielbein e4, commute with all the bosonic 
dynamical variables, especially with p's that are canonical momenta of the coordinates z's. The 
remaining canonical commutation relations are 

{^t^} = 5 AB 5 aP . (16) 
The Lagrangian (^) is invariant under the N=4 supersymmetry transformations, 

6 (0) z* = er, (17) 

5 (0) r = ~iz^°e - - l5 G»e, (18) 

d (a) z»=l^ u e {a) r, (19) 

6 {a) (I^ „r) = -«V Ha) ~ K x l {a > p l^ x a e (a) rr - 7 5 G^ e (a) , (20) 

where e and e( a ) are spinor parameters. In order to obtain supercharges, we define supercovariant 
momenta by 

Kv=P»-^AB^ A l°^ B , (21) 

where uj a b p is the spin connection. The corresponding N=4 SUSY generators in real form are 
then 

Qa=^ + i(l°l^)aG^ (22) 

q£0 = v ^ + *( 7 75 ;rW „^%G M , (23) 



5 



which satisfy the following SUSY algebra with a central extension: 

{Qa, Qp} = {Q<£\Q ( p a) } = 26 a p H + 2i( 7 °75W Z, (24) 

{Q a ,Q$} = 0, (25) 

{Qi a) ,Q?} = (a/ b). (26) 



The Hamiltonian TC and the central charge Z read 

n = \ (-^^y/ggr-Kv + - -iw^V'yyy + d^g^^ , (27) 



2 = G% - -(I^G^'VV 1 '- (28) 
It is easily checked that the central charge Z indeed commutes with all SUSY generators. 

For spectrum analysis, SUSY generators in complex form are more useful. Introducing ip = 
-71(^1 - #2)' an d defining Q = -^(Qi - iQi), one finds 

Q = p% + i^G^ (29) 

Qt = _ i^G! Mj (30) 

which generates the following simple algebra: 

{Q,Q t } = {Q (a) ,Q (a)t } = 2W, (31) 

{Q,Q} = {Q (a) ,Q (a) } = -{Q t ,Q t } = -{Q {a)t ,Q (a)t } = 2iZ, (32) 

{Q,Q {a) } = {QlQ {a) } = 0, (33) 

{QW,Q( 6 )} = {gW^QW} = 0, (a ^ 6) (34) 



It is easy to read off the BPS condition for quantum states that preserves half of supersym- 
metries. Depending on the sign of central charge, we find 

(Q T *Q t )|*>=0, (35) 

so that the given state may saturate the condition TC = ±Z. We can express this BPS condition 
in more geometrical fashion by transcribing the wavefunction to differential forms on the moduli 



space |]13 |. Note that 



[^,^] = -r>", (36) 
K,p*] = r<>;, (37) 

{^,<p* u } = 5Z. (38) 
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Furthermore, the wavefunction has the following general form, 

l*>=E^ n m-Mp(^ 1 ---^l°> ( 39 ) 

(p*"\0) = 0. (40) 

with a inner product defined by 

<$|<&'> = / dx^gJ2( nlX1 "^T n U-^ ( 41 ) 

The coefficients £1^...^ are completely antisymmetric and may be regarded as those of a p-form. 
In this language where we interpret and 99* as a natural cobasis dz^ and a natural basis g|jr, 
one finds that the following replacement can be made: 

icfficn d , i(p*^TTfi -5, (42) 
vfGn^ia, iZ ^C G = di G + i G d, (43) 

where denotes the natural contraction of the vector field G with a differential form. The BPS 
equation now becomes 

(d-i G )n = Ti(S-GA)n, (44) 

where we use the same symbol G for both the Killing vector and the 1-form obtained by contraction 
with metric. Solving this first order system, we should recover all 1/2 BPS and 1/4 BPS states of 
the underlying Yang- Mills field theory. 



4 Angular Momentum and Eigenstates on S 3 
4.1 Supermultiplet Structure of 1/4 BPS Dyons 

1/4 BPS dyons have been constructed in several difference guises. The first was as three-pronged 
strings ending on D3 branes O, while the field theoretical construction was as exact classical 
solitons 0, Q]. Neither of these was convenient for finding their supermultiplet structures; there 



are subtleties of the respective moduli space dynamics [14, 15 1. The third method, also present 
in Ref. Q, was a perturbative one. In this setup, one assumes very small electric coupling and 
work in vacua where the 1/4 BPS dyons would be stable. The construction proceeds by finding the 
lowest quantum excitation modes around purely magnetic background. The lowest are massless 
moduli. The next lowest is massive, and turns out to induce quantized electric charges to the 
system when excited, and produces many degenerate states with the same electromagnetic charges. 
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In the simplest case of a + /3 magnetic charge in SU (3), the total degeneracy for dyons of relative 
charge q ^ is 

2 6 x \2q\. (45) 

Note that q is quantized in half-integers. The highest spin of this supermultiplet is \q\ + 1. The 
relative charge defined through the following expression of electric charge one may excite on the 
system, 

q = (n + q) a + (n — q) (5. (46) 

where 2n is an integer. The consistency with Dirac quantization condition, along with the spectrum 
of the original field theory, actually demands that n ± q are integers. Thus quantization of n is 
correlated with that of the relative charge q: a half-integral q comes with a half-integral n, and an 
integral q comes with an integral n. 

When we reconstruct the dyons as bound states in the low energy dynamics, this correlation 
naturally emerges from the form of total moduli space 

M=R*x^-^, (47) 

where the quotient action of Z is crucial. In this note, we will not dwell on this point. It suffices 
to say that all dyonic states of relative charge q can be built, provided that the relative part of the 
wave function is found. 

Once we consider the wavefunction to be a tensor product of two parts, one over the relative 
moduli space, and the other over the center-of-mass moduli space, the degeneracy is more naturally 
organized as, 

2 4 x ( (2\q\ + 1) + (2\q\) + (2|g|) + (2\q\ - 1)) . (48) 



The common factor 2 4 follows from the low energy dynamics trivially. Among the low energy degree 
of freedom, there are four bosonic and eight (real) fermionic coordinates that are associated with 
the center-of-mass motion and thus free. These eight fermions acts as four pairs of massless har- 
monic oscillators, whose excitations lead to 2 4 = 16 degeneracy and the subsequent supermultiplet 
structure of N = 4 vector multiplet. 

The rest of the degeneracy factors must arise from the relative part of the dynamics, as we 
construct the 1/4 BPS dyons as bound states of monopoles. As the above decomposition suggests, 
the spin content of the supermultiplet found in Ref. [§] is such that the bound state wavefunctions 
over the relative moduli space are in four multiplets of angular momenta; one with I = \q\, two with 
I = \q\ — 1/2, the last with I = \q\ — 1. Thus, finding these BPS bound states explicitly presupposes 
detailed understanding of angular momentum in the low energy dynamics. 
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4.2 Isometries and Symmetries 



Let us recall the geometry of the Taub-NUT manifold. The metric is 

1 + I) [dr 2 + rV 2 + rVf] + a 2 , (49) 

where the <r a 's are 1-form frame on 5 3 as in Section 2. The Taub-NUT manifold has four Killing 
vectors, three of which generate SU(2) rotation of S 3 . These SU{2) Killing vectors, which we 
denote by L a , are 16], 



L\ = — sin d>dff — cot cos dd^-l —r- <9y, (50) 

sm 6 

L2 = + cos 4>d$ — cot 6* sin (/> <9<a + — — — cL, (51) 

sm 

= 84,. (52) 
The <T a 's are easily seen to be invariant under the action of these vector fields: 

C La (a b ) = 0. (53) 
The operator J a = —i£L a satisfies the usual SU(2) algebra. 

The fourth Killing vector corresponding to internal U(l) gauge rotations of monopoles, to be 
denoted as K, is 

K = 8 X . (54) 

This is precisely the triholomorphic Killing vector field that enters the low energy dynamics. This 
vector field do rotate o\ and 02 among themselves, 

- %Lk (01 ± ia-i) = ± (a\ ± ia 2 ), (55) 

but the Taub-NUT metric itself is invariant under such rotations. Since the range of x 1S 47r, the 
eigenvalue of —iCk is quantized at half- integers. Its origin as a gauge rotation generator also tells 
us that its eigenvalue should be identified with the relative charge q Q. 

These isometries of the Taub-NUT manifold naturally generate symmetries of the low en- 
ergy dynamics. In the absence of potential terms due to G = aK, this is especially clear since 
the Lagrangian is completely determined by the metric alone. Furthermore, the geometrical in- 
terpretation of the wavefunction \Q) as differential form suggests that these symmetries act on 
f2 geometrically; In fact, it is easy to see that the Noether charge associated with each of these 
isometries do acts on wavefunctions/differential forms as Lie derivative. Thus, for the rotational 
SU(2), the above operator J a may be regarded as the symmetry generators. And so is —iCk for 
the relative U(l) gauge symmetry. 
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Symmetries generated by these vectors remain symmetries of the low energy dynamics when 
we turn on the potential term determined by G = aK. This is guaranteed because L a 's and K 
preserves G: 

J a (G) = = £ K (G). (56) 

However, the J a 's are not quite the physical angular momentum, it turns out. We will come back 
to this crucial point shortly 



4.3 Angular Momentum on S 3 and Spin 

As we observed from the known degeneracy of 1/4 BPS dyons, the bound state with relative charge 
q must have angular momentum that scales with q linearly. How do we realized such a multiplet? 
Note that the Taub-NUT space has the topology of i? 4 = R + x S 3 . On the S 3 , the SU(2), which 
itself is topologically S 3 , acts freely as translations, so functions on S 3 naturally fall under various 
SU (2) representations. The angular momentum eigenfunctions on S 3 are well-known to those 
familiar with classic angular momentum theory, and are often denoted by D J mk [17]. Abstractly it 



is defined as a finite rotation operator sandwiched between a pair of eigenstates of total angular 
momentum j, 

VLM>X) = V2J+1 ti',k\U(0,<l>,x)\j;m), (57) 

where we chose a normalization so that the norm is independent of j, m or k. The definition makes 
it clear that m and k are bounded below and above by ±j, 

j >m,k> -j, (58) 

and separated from j in integer steps. For any fixed k or m, the 2j + 1 functions indexed by m or 
by k, would form a spin j multiplet of SU(2)l or SU(2)r, respectively, where the isometry group 
of undeformed S 3 is 50(4) = SU(2) L x SU(2) R . 

One might worry that the angular momentum eigenstates of our low energy dynamics may 
have little to do with these D-functions. After all, a wavefunction should be regarded as a differential 
form whose transformation properties are generically more complicated than functions. However, 
as we pointed out earlier, the Taub-NUT manifold admits orthonormal basis w^ 1 , all four of which 
are invariant under J a . This means that as long as we construct the wavefunction/differential form 
in the basis spanned by the orthonormal frame uj^, its transformation property under J a originates 
entirely from its coefficient functions. Thus we expect the generic form of the wavefunction could 
be written as, 

n = E vUP* & X) Ai;fc(r; dr, a a ), (59) 
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where the differential forms Aj-k have no explicit dependence on the three Euler angles except 
through the cr a 's. D J mk with m = —j, + form a multiplet under J a with — id x charge k. 

Given the definition of D functions, it is not difficult to show that the other lower indices, 
namely k, is an eigenvalue of the operator —id x . (Because the actual geometry of the three sphere 
is deformed, the SU (2) under which D J mk , k = —j, + would have formed a representation 

is no longer a symmetry.) Since j is bounded below by \k\, a large relative charge necessarily implies 
a large angular momentum; The degeneracy has to scale linearly with increasing k ~ q. This is 
precisely the behavior we saw from the state counting of subsection 4.1. 

However, there is something missing. Given a single eigenstate f2, we expect to generate other 
physical states related to it by acting with a SUSY charge, say Q. The underlying field theory 
tells us the new state must be fermionic/bosonic if |0) is bosonic/fermionic. In particular 

their physical spin should differ by 1/2. On the other hand, Q itself is invariant under the action 
of J a = —iCi a , and cannot impart additional angular momentum quantum numbers, it seems. 

The resolution of this dilemma is that we should modify the angular momentum operator 
by adding a "spin" piece. The hyperKahler structure of the Taub-NUT manifold supplies such 
additional conserved quantities, fortunately |1£]. Let be the three complex structures as 



before. Define a triplet of operators S a acting on fermions by, 



S a = |Z<SW- (60) 



The pointwise action is, in geometrical terms, 



S a {dz») = -T {a) \dz\ (61) 

s * il) - {w) ■ 

which is nothing but the action of the three complex structures up to a numerical factor. These 
S a 's span an SU(2) R-symmetry of the N = 4 superalgebra. 

The S a 's themselves form a triple under the JaS, and using this fact we may write down a 
new set of angular momentum generators, 

M a = J a - S a , (63) 
which also span an SU(2) algebra. Note that, unlike J a , M a commute with the R-charges S a . 

We seem to have two possible choices of angular momentum; J a which rotate R-charge and 
M a which do not. Both commute with the Hamiltonian, but their commutators with supercharges 
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are quite different matter. Under J a , the four complex supercharges fall into a singlet plus a triplet, 
i.e. 

[Ja, Q] = 0, [Ja, Q(b)} = ^abcQ(c) ■ (64) 

On the other hand, since the four complex supercharges belong to doublets under S a , they must 
form doublets under M a as well. More specifically, the following linear combinations 

Q + = Q + iQ (3) , 

Q-=iQ {1) +Q {2) , (65) 
form one doublet under M a with [S3, Q±] = — \Q±, and the second combination 

Q-=Q-iQ {3 \ (66) 

form another doublet under M a with [S3,Q±] = +\Q±- Since supercharges should carry phys- 
ical spin 1/2, we surmise that M a rather than J a should be interpreted as the physical angular 
momentum. We will denote the eigenvalues of M 2 by 1(1 + 1). 



4.4 Anti-Self-Dual Ansatz 



Now we may proceed to write down the ansatz for dyonic BPS bound states. The BPS equation 
is easily seen to be invariant under the Hodge dual operation on the wavefunction. This property 
can be used to separate the self-dual part from the anti-self-dual part of the trial wavefunction, so 
a BPS wavefunction should be either self-dual or anti-self-dual. 

Does the dynamics prefer one to the other? Hamiltonian has three kinds of potential terms. 
In addition to the purely bosonic potential G 2 , there are two more terms; one is a fermion bilinear 
contracted with dG, while the other is a fermion quadrilinear contracted with the Riemann cur- 
vature. The salient point is that dG and the Riemann curvature are both anti-self-dual tensors 
on Taub-NUT manifold. Because of this, a self-dual ansatz will not be sensitive to some spin-spin 
type long range interaction, which could be crucial for the formation of bound states. In fact, the 
threshold bound state of SU(3) monopoles (when G = 0) is known to be anti-self-dual, while no 
such self-dual bound state exists |§ . We expect that this behavior persists when G = aK is turned 
on, which motivates us to look for 1/4 BPS dyonic bound states with an anti-self-dual ansatzQ 

Let us start with anti-self-dual 2-forms. One interesting property of anti-self-dual 2-forms on 
hyperKahler 4-manifolds, is that they are of type (1, 1), upon Hodge decompositions with respect 
4 Nonetheless, there is no reason to preclude self-dual bound states that do not saturate the BPS bound. 
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to any one of three complex structures [19|. Recall that, up to a numerical factor, the "spin" S a 
act on differential forms as complex structures do. Since any form of type (n, n) is annihilated by 
the complex structure, we conclude that the 2-form part of an anti-self-dual ansatz carries no spin. 

On the other hand, the BPS equation connects even forms with even forms, so an ansatz 
containing anti-self-dual 2-form may contain, in addition, 0-form and 4-form. Neither carries "spin" : 
0-form is obviously invariant under S a , and a 4-form is always proportional to the volume form 
which is always of type (n,n). Thus, an anti-self-dual even form is always "spinless"; I equals j. 
The angular dependence of anti-self-dual even forms can be written entirely in terms of D functions 
and of the basis <r a 's. 

Of the four angular momentum multiplets, the cases of I = q > and of I = q — 1 > belong 
to this category. For I = q, the wavefunction should have the form, 

^m;q = ^mq^q;q + -^^(g-l)^ 1 ?;?-!' (6?) 

for any value of m = q, q — 1, . . . , —q, where differential forms A's can be written entirely with r 
and uj^ only. This state has relative charge q > if and only if the anti-self-dual form A's satisfy 

- i£ic\;q = 0, -iC K ^q;q-l = A^q-l- (68) 

Since the only charged combination one can build out of r and is uj 1 ± iu> 2 , which has ±1 charge 
respectively, the A's may contribute ±1 to the total charge q at most. Furthermore, only u + iuj 2 
may enter since D^,^ , which should accompany co — iu> 2 , does not exist for positive q. 

These considerations constrain possible form of A's quite severely. We find that only the 
following choice is consistent with the known quantum numbers of 1/4 BPS dyons, 

A 9;9 = f(r) + h(r)(u° A^+u 1 Aw 2 )+/(r)(w° Aoj 1 Au 2 Au 3 ), (69) 
A 9;3 _i = (b(r)/r)(u; +iuj 3 ) A^+ico 2 ), (70) 

for I = q > BPS state.0 

In case of I = q — 1 with q > 1, the anti-self-dual ansatz is even more restrictive. The only 
such ansatz consistent with known quantum numbers is, 

= p(r) D^ q _ x) (u, + ^ 3 ) A (u 1 + ^ 2 ), (71) 

for m = q — 1, q — 2, . . . , 1 — q. No 0-form may appear since a factor of (w 1 + ito 2 ) is necessary to 
make the electric charge to be q, while 4-form is ruled out subsequently by the fact that the BPS 
wavefunctions we are looking for are all anti-self-dual. 

5 We have defined the Hodge dual operation with respect to the volume form, — w° Alu 1 /\uj' 2 A cj 3 . 
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We will postpone discussion of the two remaining multiplets of I = q — 1/2 > to the following 
section. The corresponding multiplets are in odd forms, which can be found by acting supercharges 
on and ^T^. For these two multiplets, S a contribution to the physical angular momentum 
does not vanish. 



5 Dyonic BPS Bound States 

Here, we will solve for the dyonic BPS bound states explicitly. Such dyonic bound state do not 
exist for all relative charge q. Rather, it is known that \q\ must be smaller than the critical charge 
Igcrl [0, |, Hi, where 

q cr = lim a(K, K). (72) 

r— >oo 

With our current normalization, ', asymptotes to 1 at infinity, so the critical charge q cr is 
equal to the parameter a. Thus we expect to find the BPS bound state, that is, a normalizable 
and regular wavefunction that preserves half of low energy supersymmetry, only when \q\ < \a\. 
Without loss of generality, we will take both a and q to be nonnegative. 

Classical analysis of Ref. (l|, || leaves it unclear whether the bound state should exist (at 
threshold) when a = q. Classically two monopoles are infinitely separated, but quantum mechan- 
ically, there may be a bound state with powerlike decay. As the following analysis will show, 
however, no such threshold bound state exists, except for a = q = case. 

5.1 The I = q - 1 Multiplet 

The case of I = q — 1 is the simplest. Starting with the ansatz, 

n^J = p(r) (w° + *u, 3 ) A (oo 1 + iu 2 ), (73) 

the BPS equation reduces to a single ordinary differential equation for p(r), 

— r p(r) = — Air) r p(r) , (74) 
dr 

where the quantity A is defined to be the following combination: 

A = a-q(l + ^\. (75) 

A useful fact we employed is that 

e r (cj° + iu 3 ) A (a/ + ico 2 ) (76) 
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is a (nonnormalizable) harmonic 2-form. The equation is easily solved to give the wavefunction: 

n£3 = D M q -xf q ' le ' {a ' q)r + A + ^ 2 )- ( 77 ) 

The wavefunction is exponentially small at large r and normalizable as long as a > q. This way, 
we have recovered the fact that a is the critical electric charge, beyond which no bound state may 
exist. Furthermore, the wavefunction is regular as long as q > 1, which is also consistent with the 
fact that I = q — 1 BPS bound state exists only for q > 1. 



5.2 The I = q Multiplet 



The ansatz for I = q > is a bit more involved; 



n 



q;q 



D q mq (f(r) + h(r) (uj° A u 3 + w 1 A w 2 ) + f(r) (uP A cj 3 A cj 1 A cj 2 ) 



+ (Wr)/r) («° + io; 3 ) A (a, 1 + iu, 2 ) 

under which BPS equations reduce to 



d , b d , 2/t 

dr r z dr 1 + r 



•U-3, I' '/(/-/'). 



dr 

In order to solve (|79| ) for general g > 0, we proceed as follows. By substituting, 
f = u (r) e -J A , h = v{r)e-J A , b = w(r)^ e -J A , 



to (|79|), one obtains 
d 



(it — w) = 0, 



d 1 



d 



(1 + r ) 2 v = 0, 



d 



(w/A) = u — v . 



(78) 



(79) 



(80) 



(81) 



dr v dr (1 + r) 2 dr dr ' 

We solve the first equation by u = Ci+w with an integration constant C\. The remaining equations 
can be combined into a single second order equation, 

' (1 + r)w 



d_ 

dr 



d / (1 + r)w 
dr" V A 



2A 



A 



2dr 



which is integrated with a second integration constant C2 to 

' (1 + r)w 



d /(l + r)u> 
dr V ^4 



2A 



.4 



+2Cir-2C 2 - 



(82) 



(83) 



Integrating this equation gives us three-parameter family of solutions. But, fortunately, the simplest 
possible solution, 

(1 + r)w 



A 



(84) 
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turns out to be the only regular and normalizable solution. Using it to generate other radial 
functions, we find the I = q bound states: 



lc, m;q 



D q 

mq 



D 



r q e ~{a-q)r 

l + r 



a + (a + 



1 



)(o>° Auj 3 +UJ 1 Aw 2 ) +auj° Aw 1 Aw 2 Aw 3 ) 



re 



-(a-q)r 



l + r 

q/2 (cj° + iw 3 ) A (w 1 +iu 2 ). 



'm(g-l) jTjT^ Ttw/zHu, -rtu/y. (85) 

Again, we find that the wavefunction is normalizable as long as a > q. The solution is regular at 
origin for all nonnegative q. We have recovered the I = q multiplet of 1/4 BPS dyon of charge 
q > 0. 



The case of q = is a bit special, where the BPS state is a purely magnetic bound state of 
the two monopoles. For q = 0, the second D function does not exist, and we must solve a modified 
BPS equation. Nevertheless, the actual wavefunction is also obtained by taking q = limit of the 
above result: 



O 

"0:0 



l + r 



a + (a + 



1 



l+r 



) (u)° Aw 3 + w 1 Aw 2 )+aa)°Aw 1 Aw 2 A w 3 ) 



(86) 



This BPS state actually preserves all supercharges of low energy dynamics. In fact, this is the 
lowest lying state of this low energy effective theory. In the limit of aligned vacua (a = 0), this 
state also reverts to the threshold bound state of two monopoles found in Ref. ||, as it should. 

5.3 The I = q - 1/2 Multiplets 

The remaining two multiplets of I = q — 1/2 can be found most easily by acting Q on Sl m . q and 
fi^Tq found above. From O^J, we find 2q — 1 states 



,<j -(a-q)r 



r y e 



s/r + r 2 

while £l m . q produces 2q + 1 states 

r q e -(a-q)r 



(, 1 + lW 2 )A(l +W °A W 3 )D5;_ 1) . 



(87) 



Vt + r 2 



(cj°+^ 3 ) A (1+w 1 A cu 2 )^!)^ + i{uj l +ioj 2 ) A (l+w° A w 3 )^^ 

These account for all 4q states in the two I = q— 1/2 multiplets. (Alternatively we could have used 
instead of Q. Since Q^.^ and Q^.g are invariant under complex structures, and since are 
essentially Q rotated by the complex structures, the resulting 4q states will be simply the above 4g 
states with complex structure I a acting on them.) 

However, because Q is a singlet under J a , the above states form multiplets under J a instead 
of the physical angular momentum M a . To reconstruct M a multiplets, we recall the relationship 



J a = M a + S a 



(89) 
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Since the two SU (2) generators, M a and S a , commute with each other, J a multiplets are constructed 
from M multiplets and S multiplets by the rule of angular momentum addition. On the other hand, 
we actually need to reconstruct M a eigenstates from J a and S a eigenstates, for which we need to 
reverse the procedure. Without delving into details of the computation, we present the two I = q — ^ 
multiplets. Because S3 commutes with M a , we can label the two multiplets by its eigenvalue S3. 
The first has S3 = 1/2; 



(«(+)) 



(94) r e 



1p-(a-q)r 



m;q 



1 , ..2,. h + m + k D „ 



2q + l 



- m 



2q(2q + l) 



2nH 

2q (™4)(<K) 



The other multiplet has S3 = —1/2; 



(«(-)) 



(9-s) 



r Q e -(a-q)r r 



y/r + r 2 



(uP+iu?) A (lW A u?)J q n m+ ) Df ^ 
v ; v '\ 2q + l ( m -\)i 



+i{u 1 + iu: 2 ) A (1+,° A lj 3 )(J } ,_ m + i Dl 



' q + m 



1 



2q(2q + 1) V 
In both expressions, the index m takes values q — 1/2, q — 3/2, ...,—<? + 1/2. 



2 n«-l 

2g («^i)(^i) 



(90) 



(91) 



In a direct construction of these two / = q — | multiplets from the Z = g or I = q — 1 
multiplets, the role of the doublet supercharges in ( p5| ) and fl66| ) can be easily identified. As a 
simple application of the angular momentum addition rule, the operations of the doublets on the 



, 



or (g-|)©(g-|) =i®(«-i) 



+ I and I = q — | multiplets vanish identically. Since the 



even- form multiplets will produce (g + |) © (g 
but one may check that both I = q + i 
doublet Q± carries the spin eigenvalue S3 = —5, the operation of this doublet on Z = g or Z = q — 1 
multiplets, produces the I = q — \ multiplet with S3 = — \. Similarly, the application of Q± results 
in the Z = q — 5 multiplet with S3 = +i. 



5.4 Characteristics of the BPS states 



In the construction of the BPS bound states, we have limited ourselves to the case of nonnegative 
electric charge. For negatively charged bound states, we note the fact that the complex conjugation 



of a solution to the BPS equations in (44), gives another solution. Both the eigenvalues of charge 
and M3 reverse their signs under the complex conjugation. Thus the negatively charged solution, 



(q > 0) is simply given by the complex conjugation of Q l 



-m — q 



m;q- 



We now turn to the case that the Higgs misalignment parameter a is negative. When a is 
replaced by —a, only the fermionic term that couples to the Killing potential changes its sign in 
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the supersymmetric Lagrangian (|l2|). By the parity operation — > or equivalently ip — ► 
one can bring this Lagrangian to the original form with a > 0. The corresponding transformation 
amounts to replacement of dx^ by idx^. The solutions for —a are obtained if one replaces all lo^ 
of the above BPS solutions by iuA These exhaust all the possibilities. 

For the remainder of the section, we would like to comment briefly on some other aspects of 
the BPS states. With an excited charge q, the effective potential at large relative separation tends 
to (a 2 + q 2 )/2. Since the energy eigenvalue of the BPS multiplets is \aq\, one finds that binding 
energy of the dyons is 

-^binding = ^ ; (92) 

which tends to zero as the charge approach its critical value, as one expects from its classical 
counterpart. Another characteristic is the separation of the two monopole cores. In the classical 
limit, the separation between the two cores are given by r eq = |<?|/(|a| — \q\)- We expect the vacuum 
expectation value, 

to approach r eq in the classical limit. For a given supermultiplet with charge q, the expectation 
values are found to be dependent upon the angular momentum quantum number /. For instance, 
we find the expectation value 

< r > = a™( 1 + i) = r -( 1 + 5R)' <94) 

for the I = \q\ — 1 multiplet, and 

< r) = a™( 1 + ^) = r -( 1 + ^)' <95) 

for the I = \q\ — 1/2 multiplets. To restore the Planck constant h, we simply observe that classical 
charge q has the same dimension as \ffi. Since the difference between (r) and r eq scales inversely 
with q or a = q cr , it has to scale linearly with \ffi: 

(r) - r e 



eq 



o{Vh) . (96) 



Thus, (r) indeed approaches r eq = |g|/(|a| — \q\) in the classical limit. 



6 Conclusion 



In the low energy dynamics of 1/2 BPS monopoles, the misaligned Higgs vacua induces an attractive 
potential between monopoles of distinct types. This potential is crucial in formation of new dyonic 
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bound states of monopoles, some of which preserve 1/4 of supersymmetries in field theoretical 
sense, or equivalently, 1/2 of supersymmetries of the low energy dynamics of monopoles. Starting 
from the full N=4 supersymmetric low energy effective Lagrangian, we expressed the BPS equation 
of the system in the language of the differential form, which was then further reduced to a set of 
coupled first-order ordinary differential equations. These equations were solved analytically, giving 
1/4 BPS dyons as quantum bound states of two distinct SU(3) monopoles. Along the course of 
the construction, we have given a full account of the supermultiplet structures of the quantum 1/4 
BPS dyons. 

In this note, we have focused on the BPS saturated states in pursuit of the 1/4 BPS dyons. 
However, it is expected that there exist spectra of other dyonic bound states that do not saturate the 
BPS bound. Problem of finding these non-BPS bound states should be quite involved. Nevertheless, 
there are some additional information that might be of help. Supersymmetric sigma models on 
the Taub-NUT geometry are known to allow additional conserved quantities of the Runge-Lenz 



type Jl_6|, |lql - It seems quite plausible that this new symmetry generalize to the present low 
energy dynamics with potential. We have checked that the purely bosonic part indeed admits such 
conserved quantities. Such additional symmetries might be useful in finding the excited non-BPS 
bound states. 

Another aspect of the dynamics we did not discuss here is the scattering of dyons. The 
supersymmetric quantum mechanics we used can be thought of low energy dynamics for 1/4 BPS 
dyons in two ways. First, it produces these dyons as bound states. Second, it provides a frame- 
work where interaction among these dyons can be studied in a quantum mechanical setting. The 
dynamics is admittedly more involved than the usual moduli space dynamics, given the presence 
of the potential. It requires further study. 
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